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MEASURING METHOD FOR DECIDING DIRECTION TO A FLICKERING 
SOURCE 

5 

TECHNICAL FIELD 

The present invention relates to a method for deciding 
the direction to a flickering source in relation to a 
measurement point in an electrical network with 
10 alternating current with a network frequency with low- 
frequency amplitude variations from a flickering source. 
The present invention also relates to an arrangement for 
carrying out the method. 

15 BACKGROUND ART 

During the production of electricity, generators are used 
that produce an alternating voltage around a certain 
frequency. The user of the produced electricity is aware 
of which frequency is involved and which voltage is 

20 provided on the electricity network. The user of the 
electricity network requires as clean electricity as 
possible. That is, the user requires electricity that is 
well-defined at the specified frequency and the specified 
voltage. In normal electricity networks that supply towns 

25 and factories, it is, however, normal for the electricity 
network to comprise low-frequency voltage distortions. 
Voltage distortions arise when various devices are 
cyclically connected into the network, which devices 
comprise, for example, capacitances and inductances. 

30 

Large power-hungry installations such as induction 
furnaces, compressors, lift motors, pumps, etc., act to 



increase the amount of voltage distortion in the 
electricity network. These loads often cause low- 
frequency (<25 Hz) and periodic fluctuations in the 
effective value of the voltage. The phenomenon is usually 
5 called flicker. The most noticeable effect of flicker is 
the annoyance that can be experienced when the light 
intensity from incandescent lamps fluctuates in line with 
the voltage variations caused by the flicker. Studies 
have shown that people are particularly sensitive to 

10 light fluctuations with a repetition frequency in the 
range 0.5 Hz to 25 Hz. At the level of maximum 
sensitivity (approx. 9 Hz), the relative voltage change 
needs only to be 0.25% for people to be able to perceive 
that the light from the incandescent lamp flickers. 

15 Problems with flicker arise primarily in areas with very 
heavy industry (ironworks and paper-mills, etc.) but can 
also arise in areas with weak electricity networks and in 
the vicinity of wind-power installations. 

20 In the event of a dirty electricity network comprising 
flicker, it is of interest to know where the flickering 
source is located. An electricity producer wants to show 
that it is a consumer that makes the electricity network 
dirty and can, in such a case, demand that the consumer 

25 pay a penalty charge or corrects the problem to eliminate 
the flicker. A consumer, on the other hand, wants of 
course to show that it is not the consumer causing the 
problem. In addition, the consumer wants to have value 
for money and thus wants to show that it is the 

30 electricity producer that is supplying dirty electricity. 
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Correcting the problem of flicker is often an expensive 
operation where all or parts of the network can need to 
be reconstructed in such a way that the internal 
impedance of the electricity network is reduced, for 
5 example by use of new and larger cables. Another way to 
correct the problem is to install counter measures 
against the flicker in question. Such counter measures 
are usually very expensive. An example of a counter 
measure is "Static Var Compensation (SVC)" which 
10 dynamically controls changes in the system. 

Methods exist for determining the content of voltage 
flicker, which are described in the standard IEC-61000-4- 
15. Measuring instruments that record voltage flicker 

15 according to this standard show the presence of voltage 
flicker by calculating and displaying the parameters Ifl, 
Pst and Pit. There is, however, nothing in the standard 
or in existing measuring instruments that shows the 
direction to the flickering source in relation to the 

20 measurement point. 

There is thus a great need for a method and an 
arrangement that can determine the direction of 
propagation of flicker. The direction of propagation of 
25 the flicker shows whether the flickering source is above 
or below a measurement point and can thus be of great use 
when the flickering source is to be traced. 

DISCLOSURE OF INVENTION 
30 The present invention intends to solve the problems that 
are described above of low- frequency flicker that gives 
rise to periodic fluctuations in the effective value of 
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the voltage. In the following, such periodic fluctuations 
are called low-frequency variations or flicker. 

The problem is solved by a method for deciding the 
5 direction to a flickering source in relation to a 
measurement point in an electricity network with 
alternating current with a network frequency f c with low- 
frequency amplitude variations from the flickering 
source. The method provides that an amplitude-modulated 

10 current signal and an amplitude-modulated voltage signal 
are recorded. Both the amplitude-modulated current signal 
and the amplitude-modulated voltage signal are signal 
processed in such a way that the low- frequency amplitude 
variations in both the current and voltage are separated 

15 from the carrier wave in the form of a flicker component 
for current and a flicker component for voltage. The 
flicker component for current is thereafter multiplied by 
the flicker component for voltage in such a way that a 
product is produced. The product is processed in such a 

20 way that a flicker power IT is obtained with a sign value 
that indicates in which direction a flickering source is 
located in relation to the measurement point. 

According to an embodiment of the invention, the sign 
25 value of the flicker power at the measurement point is 
negative when the flickering source is located below the 
measurement point and positive when the flickering source 
is located above the measurement point. 

30 The method has two important advantages : 
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1. The flicker power at the measurement point is 
determined. By this, the flickering source can be traced. 

2. The method is also accurate when the measurement 
5 is carried out via voltage and current transformers 

installed in the network. The frequency spectrum of a 
signal with flicker consists of a carrier wave {for 
example 50 or 60 Hz) and of side-bands on both sides of 
the carrier wave, with frequency distance to the carrier 

10 wave corresponding to the flicker frequencies. The 
frequency spectrum of the whole signal package (carrier 
wave and the low-frequency signals) is thus in a narrow 
frequency band around the carrier wave, which means that 
measurement via existing current and voltage transformers 

15 can be used as these are designed for highest accuracy 
around the network frequency. 

The invention can be carried out by a number of methods 
described below. 

20 

In a first method (method 1) according to an embodiment 
of the invention, the method includes the steps of: 

- recording of an amplitude-modulated current signal i(n) 
comprising signals that originate from the network 

25 frequency f c and the low-frequency amplitude variations 
in the current signal i (n) ; 

- recording of an amplitude-modulated voltage signal u(n) 
comprising signals that originate from the network 
frequency f c and the low-frequency amplitude variations 

30 in the voltage signal u(n); 

- signal processing of the current signal i(n) in such a 
way that only the low-frequency amplitude variations 
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remain in the form of a flicker component for the current 
signal i (n) ; 

- signal processing of the voltage signal u(n) such a way 
that only the low-frequency amplitude variations remain 

5 in the form of a flicker component for the voltage signal 
u(n) ; 

- creation of a product by multiplication of the flicker 
component for current and the flicker component for 
voltage; 

10 - processing of the product in such a way that a flicker 
power n is obtained with a sign value that indicates in 
which direction the flickering source is located in 
relation to the measurement point. 

15 An advantage of the first method is that it does not 
require extensive computer capacity, but can easily be 
implemented in a suitable device. 

According to another embodiment of method 1, the method 
20 comprises : 

the signal processing of the current signal i(n) 
comprising the steps: 

creation of a first demodulated signal by 
demodulation of the current signal i (n) ; 
25 - filtering off (eliminating) of the signals that 

originate from the network frequency f c in the first 
demodulated signal in such a way that only the low- 
frequency variations remain in the form of the 
flicker component for current; 
30 - the signal processing of the voltage signal u(n) 
comprising the steps: 
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creation of a second demodulated signal by 
demodulation of the voltage signal; 

- filtering off (eliminating) of the signals that 
originate from the network frequency f c in the 
5 second demodulated signal in such a way that only 

the low-frequency variations remain in the form of 
the flicker component for voltage. 

According to another embodiment of method 1, the method 
10 comprises the steps: 

- filtering off of the signals that originate from the 
network frequency f c in the first demodulated signal in 
such a way that only the low- frequency . variations 
relating to the flicker component for current remain in 

15 the form of a flicker signal Ilf (n) for current; 

- filtering off of the signals that originate from the 
network frequency in the second demodulated signal in 
such a way that the low-frequency variations relating to 
the flicker component for voltage remain in the form of a 

20 flicker signal U LF(n ) for voltage; 

the product creating an instantaneous power signal 
ri(n) by multiplication of the flicker signal Ilf (n) for 
current and the flicker signal U LF(n) for voltage; 

- the product being processed by the creation of the 
25 average value of the instantaneous power signal Il(n) 

whereby the flicker power n is created with the sign 
value . 

According to yet another embodiment of method 1, the 
30 method comprises: 

- the first demodulated signal being created by square 
demodulation of the current signal; 
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- the second demodulated signal being created by square 
demodulation of the voltage signal. 

According to an embodiment of method 1, the filtering is 
5 carried out with a bandpass filter with a lower limit of 
0.1 Hz and an upper limit of 35 Hz. A preferred upper 
limit is, however, 25 Hz. 

In a second method (method 2) according to an embodiment 
10 of the invention, the method includes the steps of: 

- recording of an amplitude-modulated current signal i (n) 
comprising signals that originate from the network 
frequency f c and the low-frequency amplitude variations 
in the current signal i (n) ; 

15 - recording of an amplitude-modulated voltage signal u(n) 
comprising signals that originate from the network 
frequency f c and the low-frequency amplitude variations 
in the voltage signal u(n); 

frequency analysis of the wave form of the voltage 

20 signal u(n) by an N-point DFT analysis {Discrete Fourier 
Transform) , which gives rise to a voltage vector U which 
contains the frequency spectrum for the voltage signal 
u(n) in the form of N complex voltages; 

- frequency analysis of the wave form of the current 
25 signal i (n) by an N-point DFT analysis (Discrete Fourier 

Transform) , which gives rise to a current vector I which 
contains the frequency spectrum for the current signal 
i(n) in the form of N complex currents; 

- creation of a power vector P by the multiplication, 
30 element by element, of the voltage vector U and the 

current vector I; 
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- multiplication of the power vector P by a weighting 
vector W that eliminates the power component that 
originates from the network frequency, with the power 
vector P comprising partial powers P k concerning power 

5 components from the flickering source, 

- creation of a flicker power n with a sign value by 
summation of the partial powers P k , and 

analysis of the sign value, with the sign value 
indicating in which direction from the measurement point 
10 the flickering source is to be found. 

In the second method, the voltage signal u(n) and the 
current signal i (n) are signal processed by the frequency 
analysis described above. The product given by method 1 

15 corresponds to the creation of the power vector P in 
method 2. The flicker components for current and voltage 
in method 1 do not have anything directly corresponding 
to them in method 2, but flicker components arise in the 
form of power components P k concerning power components 

20 from the flickering source after the power vector P has 
been multiplied, element by element, by the weighting 
vector W which comprises the elements w k . The processing 
described in method 1 corresponds to the summation of the 
partial powers P k in method 2. 



According to an embodiment of method 2, the flicker power 
n is created by the following step: 

- summation of the partial powers P k by the formula: 



25 
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An advantage of this method is that no demodulation 
residues arise. 

According to another embodiment of method 2, the flicker 
5 power n is created by the following steps 

- square demodulation x 2 of the voltage signal u(n); 

- square demodulation x 2 of the current signal i (n) ; 

- calculation of the frequency spectrum of the square- 
demodulated voltage signal by an N-point DFT analysis 

10 (Discrete Fourier Transform) which gives rise to the 
voltage vector (U) ; 

- calculation of the frequency spectrum of the square- 
demodulated current signal by an N-point DFT analysis 
(Discrete Fourier Transform) which gives rise to the 

15 current vector (I); 

- creation of the flicker power n by summation of the 
partial powers P k which contribute to the flicker 
phenomenon by the formula: 



20 where the elements wl k and w2 k replace W k and eliminate 
the power component that originates from the network 
frequency and weight the correct amplitudes of the 
frequency components U k and I k , in accordance with: 
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w\ k = 



— for \<k<i 
U c 

0 fork> i 



w2 k = 



— for \<k<i 
h 

0 for k> i 



where it is assumed that the low-frequency flickers are 
to be found in a frequency band up to and including tone 

5 i (0 < f flicker < i) • 

As shown above, both method 1 and method 2 give rise to 
the flicker power n according to the invention with a 
sign value that indicates in which direction a flickering 
10 source is located in relation to a measurement point. It 
is thus possible to signal process the current signal and 
voltage signal in both the time plane and the frequency 
plane in order to obtain the required flicker power n. 

15 The invention also relates to an arrangement comprising 
for carrying out the methods described above. 

According to an embodiment of the invention, the 
arrangement comprises means for deciding the direction to 
20 a flickering source in relation to a measurement point in 
an electricity network with alternating current with a 
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network frequency f c with low- frequency amplitude 
variations from the flickering source. The arrangement 
also comp rises: 

means for recording an amplitude-modulated current 
5 signal i(n) comprising signals that originate from the 
network frequency f c and the low-frequency amplitude 
variations in the current signal i (n) ; 

means for recording an amplitude-modulated voltage 
signal u(n) comprising signals that originate from the 
10 network frequency f c and the low-frequency amplitude 
variations in the voltage signal u(n); 

- means for signal processing the current signal i (n) in 
such a way that only the low-frequency amplitude 
variations remain in the form of a flicker component for 

15 the current signal i(n); 

- means for signal processing the voltage signal u(n) in 
such a way that only the low-frequency amplitude 
variations remain in the form of a flicker component for 
the voltage signal u(n); 

20 - means for creating a product by multiplication of the 
flicker component for current and the flicker component 
for voltage; 

- means for processing the product in such a way that a 
flicker power n is obtained with a sign value that 

25 indicates in which direction the flickering source is 
located in relation to the measurement point. 



According to another embodiment of the invention: 
- the means for signal processing of the current signal 
30 (i(n) comprises: 

- means for creating a first demodulated signal by 
demodulation of the current signal i(n) ; 
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- means for filtering off the signals that originate 
from the network frequency f c in the first 
demodulated signal in such a way that only the low- 
frequency variations remain in the form of the 

5 flicker component for current; 

- the means for signal processing of the voltage signal 
u{n) comprises: 

- means for creating a second demodulated signal by 
demodulation of the voltage signal; 

10 - means for filtering off the signals that originate 

from the network frequency in the second demodulated 
signal in such a way that only the low-frequency 
variations remain in the form of the flicker 
component for voltage. 

15 

According to another embodiment of the invention, the 
arrangement comprises: 

- means for frequency analysis of the wave form of the 
voltage signal u(n) by an N-point DFT analysis (Discrete 

20 Fourier Transform) , which gives rise to a voltage vector 
U which contains the frequency spectrum for the voltage 
signal u(n) in the form of N complex voltages; 

- means for frequency analysis of the wave form of the 
current signal i (n) by an N-point DFT analysis (Discrete 

25 Fourier Transform) , which gives rise to a current vector 
I which contains the frequency spectrum for the current 
signal i (n) in the form of N complex currents; 

- means for the creation of a power vector P by the 
multiplication, element by element, of the voltage vector 

30 U and the current vector I; 

- means for the multiplication of the power vector P by a 
weighting vector W that eliminates the power component 
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that originates from the network frequency, with the 
power vector P comprising partial powers concerning 
power components from the flickering source, 

- means for the creation of a flicker power II with a 
5 sign value, by summation of the partial powers P k , and 

- means for analysis of the sign value, with the sign 
value indicating in which direction from the measurement 
point the flickering source is to be found. 

10 BRIEF DESCRIPTION OF DRAWINGS 

The invention will be described below in the form of a 
number of embodiments with reference to a number of 
figures, in which: 

15 Figure 1 shows an equivalent two-terminal network for an 
electricity network according to the invention; 

Figure 2a shows variations in the effective value U L of 
the voltage; 

20 

Figure 2b shows variations in the effective value I L of 
the current; 

Figure 3 shows a frequency spectrum for an amplitude- 
25 modulated voltage signal with only one low-frequency 
component; 

Figure 4 shows a signal flow chart for a measurement 
method according to an embodiment of the invention; 

30 

Figure 5 shows a signal flow chart for a measurement 
method according to a second embodiment of the invention; 
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Figure 6 shows a signal flow chart for a measurement 
method according to yet another embodiment of the 
invention; 

5 

Figure 7 shows an amplitude characteristic for a bandpass 
filter that is used in the embodiment described in 
association with Figure 6; 

10 Figure 8 shows an outline block diagram of an arrangement 
that can be used for the methods described in the 
different embodiments; 

Figure 9 shows schematically a network comprising a 
15 flickering source Fl, a load LI, a generator G for 
generating an alternating voltage, and 

Figure 10 shows schematically a diagram of the flicker 
power FT for a number of sampling points n. 

20 

PREFERRED EMBODIMENTS 

Both the general theory and a number of embodiments will 
be described below. The general theory is necessary for 
an understanding of the embodiments described below. In 
25 the equations that are described, a point between two 
letters in an equation indicates a multiplication of 
vectors, element by element. 

As mentioned above, the invention relates to a method for 
30 determining in which direction a flickering source is 
located in relation to a measurement point. To determine 
the direction of low-frequency variations such as 
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flicker, the wave form of the voltage and current must be 
recorded in the phase or phases where the direction is to 
be determined. Thereafter, the recorded information is 
signal processed according to any one of the embodiments 
5 that are described in the following signal flow charts. 
The result of the signal processing is a flicker power n 
with a sign value. The sign value indicates in which 
direction the flickering source is located in relation to 
a measurement point. 

10 

Figure 1 shows an equivalent two-terminal network for an 
electricity network according to previously known 
technology. The electricity network can be divided 
schematically into three parts which are commonly called 

15 generator 1, transmission line 2 and load 3. When a load 
3 is connected, a current I will flow in the circuit. 
This leads to a drop in voltage U Z t arising across the 
internal impedance Z T so that the voltage U L across the 
load drops. If the load is connected in and disconnected 

20 periodically (cyclically) , the current I will also 
increase and decrease cyclically, so that the voltage U L 
decreases and increases cyclically (provided that the 
generator voltage U G remains constant) . 

25 Figure 1 shows a measurement point 17 and a point 18 that 
marks a point above, that is upstream of, the measurement 
point 17, and a point 19 that marks a point below, that 
is downstream of, the measurement point 17 . The terms 
above and below the measurement point 17 are important, 

30 as it is necessary to indicate in which direction the 
flickering source is located in relation to the 
measurement point. 
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According to the invention, a recording is taken (at the 
measurement point) of an amplitude-modulated current 
signal i (n) comprising signals that originate from a 
5 network frequency f c and low-frequency amplitude 
variations in the current signal i (n) . The low-frequency 
amplitude variations originate from the flickering 
source. In addition, a recording is carried out of an 
amplitude-modulated voltage signal u(n) comprising 
10 signals that originate from a network frequency f c and 
the low-frequency amplitude variations in the voltage 
signal u(n) . Also here, the low-frequency amplitude 
variations originate from the flickering source. 

15 Figures 2a and 2b show outline diagrams of variations in 
the effective value U RMS of the voltage for the voltage U L 
and effective value Irms of the current for the current 
voltage I. According to the invention, the changes in the 
effective values for current and voltage Urm S and Irms 

20 relate to the instantaneous changes in the voltage and 
current that arises due to the low-frequency amplitude 
variations from the flickering source. Figures 2a and 2b 
do not show signals that originate from the network 
frequency. It is, however, known that the low-frequency 

25 amplitude variations modulate the network frequency f c . 
Any changes that take place in the current and voltage 
and that originate from the network frequency are 
negligible. The effective value changes U RMS and I RMS thus 
reflect the instantaneous changes in current and voltage 

30 that can be derived from the low-frequency variations. In 
Figures 2a and 2b, from the point of view of the 
invention, it is thus only of interest to show the 
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variations in the effective value U RMS of the voltage and 
the effective value I RMS of the current. The variations in 
the effective value U RMS of the voltage and the effective 
value Irms of the current depend, according to the method 
5 according to the invention described below, on whether 
the flickering source is located above or below the 
measurement point 17. 

Figure 2a shows the situation when the flicker power is 
10 propagated from the load towards the generator, with the 
changes in the effective values for current Irms and 
voltage Urms taking place instantaneously and in antiphase 
(180 degrees out of phase) . When the load increases, the 
current I increases, whereby the drop in voltage across 
15 U ZT increases, which means that U L reduces 
instantaneously . 

Figure 2b shows the situation when the flicker power 
propagates from the generator towards the load, with the 

20 changes in the effective values for current Irms and 
voltage Urm S taking place instantaneously and in phase. An 
increase of U G gives an instantaneous increase of the 
current I, which gives an instantaneous increase of U L , 
which thus gives a simultaneous change of current and 

25 voltage. 

In a theoretical consideration of flicker, it is 
expedient to let the current signal and voltage signal 
u(t), i(t) be described as an amplitude modulation. The 
30 mathematical expression for such signals is given by [1] 
and [2], in accordance with: 
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"(0 = \U C + YjJ mk cos(<V + fi k )Jcos(© e / + p c ) [1 ] 

K0 = ^J e +^ m jr^os(a? k l + a k )jcos(co c t+a c ) [2] 

The signals consist partly of a carrier wave U c cos (co c t) 
and I c cos(co c t) respectively, and the low-frequency 
flicker signals U mk cos (co k t) and I mk cos {co k t ) respectively. 
5 The index c relates to the carrier wave's component and 
the index m relates to the low- frequency variation's 
component. The index k relates to the index for a 
summation across values k=l to N. 

10 Figure 3 shows a frequency spectrum for an amplitude- 
modulated voltage signal with only one low-frequency 
component (single-tone modulation) . The low-frequency 
component relates here to a low-frequency tone that gives 
rise to the variations in amplitude of both current and 

15 voltage. By low-frequency tone is meant here a low- 
frequency signal. 

Figure 3 shows the frequency bands to which the low- 
frequency variation gives rise. In Figure 3, f c relates 
20 to the carrier frequency and f m to the frequency of the 
low-frequency variation. The frequency components with 
single-tone modulation are to be found at the frequencies 

fc, fc-fm and fc+fm- 

25 Figure 3 also shows the amplitude U c for the voltage of 
the carrier frequency f c and the amplitude U m /2 for the 
voltage component from the low-frequency variation f c . 
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Figure 3 shows that a frequency spectrum of the 
modulating tone forms an upper side-band f c +fm and a lower 
side-band f c -f m with half the original amplitude U m in 
5 each side-band and located at a frequency distance from 
the carrier frequency f c corresponding to the modulating 
frequency f m . Low frequencies such as the low-frequency 
flicker give rise to side-bands close to the carrier 
frequency. The higher the frequency of the flicker, the 
10 larger the frequency distance of the modulating frequency 
from the carrier frequency. 

The flicker power, n, is the power that has its origin 
in the modulating tones and can be expressed as: 

n = 4 Jj ]T U nk cos(co k t + fi k ) Y ]T cosK t + a k )\lt = { orthogonally } = 

J . [3] 

t=l 2 2 

15 

The formula [3] shows that, after multiplication and 
integration, the individual low- frequency tones in 
current and voltage create the flicker power. In order to 

20 determine this power, the amplitude and phase of the low- 
frequency signals in voltage and current must be known 
and the low-frequency signals must be able to be 
extracted from the signal package [1] and [2]. This can 
be carried out in several different ways that lead to the 

25 same result. Examples of different methods are described 
below as different embodiments of the invention. 
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Two embodiments will be described below that are based on 
the frequency spectrum of the sampled wave forms u[n] and 
i[n] being determined by carrying out an N-point DFT 
analysis (DFT = Discrete Fourier Transform) . In the 
5 practical case, an FFT analysis ( FFT = Fast Fourier 
Transform) can be used, which is a calculation algorithm 
that provides the same information as the Fourier 
Transform. The output data from the analysis is two 
complex-value voltage and current vectors U and I, which 

10 contain frequency spectrums for u[n] and i [n] in the form 
of k = N complex voltages U k and currents I k . The 
frequency resolution that is obtained depends upon the 
selected sampling frequency (f s ) and the number of 
samples N that are included in the calculation, in 

15 accordance with the relationship Af = f s /N. 

Mathematically, the vectors U and I are expressed: 

«y = [t/ t , c/ 2 , t/ 3 , . . .t/ /V J = Ijc/, a , |c/ 2 1 ^ 2 , |ey 3 k A , - , |^ J [4] 

/ = [/ 1 J 2 J3,...,/ JV ,] = |/ 1 |Za 1 ,|/ 2 |Za 2 ,|/ 3 |Za3,...,|/ A ,|Z^J [5] 

The complex power, S, is created on the basis of the 

20 vectors U and I and contains the power quantities P 

(average power) and Q (reactive power) per tone in 
accordance with: 
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s = [s, , s 2 ,s 3 , ...5 iV ,] = |j£/, • /; , jU 2 ■ i; , i t/, • /; .... , | c/. v • /; j = 
-[I|y 1 |^.|/ I |z-a I> i^|^A-|/ a k-«..^,kA-|/ I k-^,.4 

-[I|Cf 1 |.|/ 1 ^.i|V 1 |.|/ > ^,I|E/,|-^ r .^ 

=[/>, +ya,p, +ye 2 ,p 3 + jQ i ,...,p N +jq. v ] 



[6] 



The active power Pi is obtained as the real part of Si, 
in accordance with: 



5 The vectors U, I and S contain the frequency spectrum 
from the angular frequency f = 0 Hz to f = f s Hz . There 
is often a desire to weight the information and/or zero 
certain frequencies. A simple and effective way of doing 
this is to introduce the weighting vector W, in 
10 accordance with: 



' The weighing factor W contains elements, w k , that contain 
constants that are multiplied, element by element, by U, 
15 I and S. If the voltage vector U is to be zeroed for 
frequencies above half the sampling frequency, the 
following operation is carried out (the point notation in 
[9] relates to multiplication element by element) : 



/»=Refo}=Rete+./e,} 



[7] 



W = [w t , w 2 ,w,,...," 




[8] 
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where 



[9] 




1 for \^k<i- i 
0 for —<k<N\ 



Umod is the modified voltage vector containing the 
5 frequency spectrum up to half the sampling frequency. 
Above this, the elements of the voltage vector are zero. 

The weighting vector's elements, w k , can be selected so 
that the required result is obtained. For example, a 
10 filter characteristic can be achieved by selecting 
suitable values of the elements w k . In addition, 
different weighting vectors can be introduced for U, I 
and S with the aim of achieving the required result. 

15 The flicker power, n, is obtained by summation of the 
powers, P k , which contribute to the flicker phenomenon. 
This is to say: 
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[10] 
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Figure 4 shows a signal flow chart for a measurement 
method according to an embodiment of the invention for a 
first measurement method. The calculation of the flicker 
power EL is carried out by frequency analysing the wave 
5 form of the voltage and current. 

The input signals for the signal flow chart in Figure 4 
consist of the two input signal vectors, u[n] and i[n], 
which contain the sampled wave forms for voltage and 
10 current. The distance in time between two indices (n and 
n+1) in the input signal vectors corresponds to l/f s , 
where f 3 relates to the sampling frequency. 

By utilizing the previously described theory for the 

15 signal flow chart that is described in Figure 4, the 
flicker power can be calculated. Figure 4 shows two 
blocks designated DFT that refer to a transformation of 
the current signal and voltage signal, i(n) and u(n) 
respectively, by a so-called Discrete Fourier Transform. 

20 The signal that comes out from the DFT for the voltage 
signal u(n) is designated U and comprises the components 
Uk^Pk which relate to a number of values for different 
angles p k / where the index k indicates which point is 
referred to in the DFT. The signal that comes out from 

25 the DFT for the current signal i (n) is designated I and 
comprises the components I k Z.q> k which relate to a number 
of values for different angles q> k , where the index k 
indicates which point is referred to in the DFT. Figure 4 
shows a multiplier 3 in which the vectors U, I and W are- 

30 multiplied. The output signal from the multiplier is 
designated P and relates to a power vector comprising 
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components partial powers P k that describe the active 
power in accordance with [7] as the real part of the 
complex power S in accordance with [6] and [7] . 

P k = Re{0.5w k - 1 U k | ZpY j U Z-<p k } 

5 [10a] 

Figure 4 shows, in addition, a summation point X that 
indicates the creation of the flicker power n by 
summation of the active powers P k in accordance with 
10 [10a], where: 

n=£ Re{0.5 w k |u k |zpk il k |z-(p k } 

[10b] 

The frequency spectrums of the wave forms u[n] and i [n] 
are calculated by an N-point DFT analysis. Thereafter, 

15 the flicker power n is calculated by utilizing [10] . The 
content of the weighting vector W is decisive for the 
final result. Elements in W are to be selected so that 
the power in the carrier wave is zeroed and so that only 
the tones that originate from flicker are included. This 

20 can be carried out in three ways. The information is used 
from either 1) the upper side-band, 2) the lower side- 
band or 3) there is a summation per frequency of the 
information in both side-bands, and the power is 
thereafter calculated by utilizing [10]. 

25 

In order to exemplify how the calculations are carried 
out, we assume that the carrier wave signal is found in 
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element M of the vectors U and I . We also assume that the 
low-frequency flicker tones in the vectors U and I are 
found in the elements k with the index: 



5 



M-i<k<M and M < k < M + i 



The number of flicker tones is given by the constant i. 
If the choice is made to utilize the information in the 
lower side-band for the calculation of the flicker power, 
10 the elements w k in the weighting vector W are selected, 
in accordance with 



The flicker power is thereafter created by utilizing the 
15 formula [10] . 

If the choice is made to utilize the information in the 
upper side-band for the calculation of the flicker power, 
the elements w k in the weighting vector W are to be 
20 selected, in accordance with 



The flicker power is thereafter created by utilizing the 
formula [10] . 
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If the choice is made to utilize the information in both 
side-bands for the calculation of the flicker power, the 
elements w k in the weighting vector W are selected, in 
accordance with 



0 for \<k<M-i 

1 for M-i<k<M and M<k<M + i 
0 for k = M 



The flicker power n is thereafter created by 
formula 



-«)•} 



10 [10c] 

Figure 5 shows a signal flow chart for a second 
measurement method according to another embodiment of the 
invention, where the calculation of the flicker power is 
15 carried out by square demodulation of voltage and 
current, after which a DFT analysis of the low-frequency 
signals is carried out. 

The method is reminiscent of the first measurement method 
20 described in association with Figure 4, with the 
difference that, before the frequency spectrum is 
calculated, the' signals are square demodulated, which is 
designated X 2 in Figure 5. By this means, the low- 
frequency signals are separated from the carrier wave, in 
25 accordance with the expression below: 
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u\t) = \u c +|x,cos(G> A V + /?,)j cos 2 + = 

= |f// +2-U c - j^U nk cos(ay + /3 k ) + cos(a> k t + fi k )J j . A (l + cos(2<y c r + 2/? e )) = 

5 r , 

C/ " 1 , ft, \2 

2 'V*'JP* cos(^+A)+^Z C/ m * cos ^+A)J |-|cos(2«t/+2/r c ) 

[11] 

Square demodulation means that two signal packages are 
created, separated in frequency. One signal package 
consists of a direct component, the modulating 
5 frequencies and the modulating signals squared. The 
latter is an unwanted combination product. The second 
signal package contains the same terms but centred in 
frequency around double the carrier wave frequency. For 
calculation of the flicker power, only the terms are 

10 included that are marked by double underlining in [11] . 
One term contains the low-frequency flicker tones 
multiplied by the carrier wave amplitude and the other 
term is the modulated signal squared. In terms of size, 
the first term is much larger than the second, which 

15 means that the squared term affects the result very 
little. Corresponding expressions are obtained when 
squaring the current signal i (t) . 

The input signals for the signal flow chart in Figure 5 
20 consist of the two input signal vectors, u[n] and i[n], 
which contain the sampled wave forms for voltage and 
current. The distance in time between two indices (n and 
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n+1) in the input signal vectors corresponds to l/f s . In 
X 2 each element in the input signal vectors is squared 
and constitutes input data for the N-point Discrete 
Fourier Transform, DFT. The output signal from the 
5 respective DFT is complex-value vectors U and I 
containing frequency spectrums for u[n] and i [n] (both 
amplitude information and phase information) with a 
frequency resolution of Af = f s /N. The content in U and I 
is expressed by: 

u =[^,^ 2 ,c/ 3 ,...^J=|j^JzA,|c/ 2 |z/? 2 ,|f/3[z/?3,...,|f/ / ,|z^,] 

/ = [/ l5 / 2 ,/3,...,/ / ,,]=|/ 1 |Z« l? |/ 2 |Z« 2> |/ 3 |Z«3 J ...,|/ A ,|Z^J 

10 

The complex power, S, is calculated by : 

[12] 

15 The sought flicker power II is then obtained by )/ the 
formula 

[13] 

20 The elements wl k and w2 k zero the frequency components 
that do not cause flicker and weight the correct 
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amplitudes of the frequency components Uk and I k , in 
accordance with 



wl k =\U c 




— for \<k<i 



0 for k>i 



[14] 



w2 k = 



0 



for k > i 



for \<k<i 



5 



[15] 



In [13], [14] and [15], it is assumed that flicker tones 
are found in a frequency band up to and including tone i 
(0 < f flicker ^ i) • The correct value for the weighting 
10 factors can be identified by studying [11] . The method is 
not limited to the weighting factors that are given 
above, other constants in the weighting vector being able 
to be used to obtain the required filter effect. 

15 Figure 6 shows a signal flow chart according to yet 
another embodiment of the invention for a third 
measurement method, where the flicker power is calculated 
by square demodulation of voltage u(n) and current i (n) 
where the low-frequency tones are filtered by the use of 

20 bandpass filters 2A and 2B respectively. 

Instead of using DFT analysis and weighting vectors as 
described in the two previous embodiments described in 
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association with Figures 4 and 5, in this measurement 
method bandpass filters 2A, 2B are used for filtering the 
low-frequency flicker tones. 

5 The input signal vectors u[n] and i[n] contain the 
sampled wave forms for voltage and current. In the 
components 1A and IB, the carrier wave (signal from the 
network frequency) and the low-frequency flicker tones 
are separated, similarly to the measurement method 

10 described in association with Figure 5, by each sample 
being squared. Only the low-frequency flicker tones are 
allowed to pass through the bandpass filters 2A and 2B. 
As input signals to the multiplier 4, there are thus only 
the low-frequency flicker tones in voltage u LF [n] and 

15 current i LF [n]. The output signal from the multiplier 4 
is the instantaneous flicker power FI[n] = u LF [n] x i LF [n]. 
The flicker power Tl is obtained by the integrator 5 
creating an average value for the instantaneous power 
ll(n). This can be carried out with a digital filter in 

20 the form of, for example, a low-pass filter. 

The bandpass filters 2A and 2B are dimensioned to obtain 
a lower limit frequency of 0.1 Hz and an upper limit 
frequency of 25 Hz. Alternatively, the bandpass filters 
25 can be used that are defined in the standard IEC 61000-4- 
15 that describes a flicker algorithm. 

The transmission function for the filter in the IEC 
standard is expressed as 
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Mi \ - ^&>i-y \+s/co 2 

s 2 +2Zs + co 2 ' (s + l/a> 3 )(s + \/eo 4 ) 

[16] 

The coefficients in [16] are to have values in accordance 
5 with the following table. 



k = 1.74802 


X = 271-4 . 05981 


a>i=27i-9. 15494 


0)3=271-1.22535 


co 2 =27t-2 .27979 


co 4 =2ti-21 . 9 



Table 1. 

10 Figure 7 shows the amplitude characteristic for the 
bandpass filter with transmission function in accordance 
with [16] and coefficients in accordance with Table 1. 

It is also possible to select bandpass filters with a 
15 different characteristic to that described in [16] . For 
example, an Mth order Butterworth or Chebishev filter can 
be selected. 

Figure 8 shows a block diagram of the hardware for the 
20 instrument. The instrument is constructed around a signal 
processor 7 which administers the measurement information 
and carries out necessary calculations according to the 
measurement methods that are described in the Figures 4-6 
and 9. The signal processor 7 also controls the sampling 
25 process for the conversion of analogue signals to digital 
signals (the A/D conversion) . The recorded signals; that 
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is, the wave form of current and voltage, are obtained 
either from current transformers and voltage transformers 
in the network or from measurement sensors that are to be 
found in the instrument. The signal processor can be a 
5 computer or a logic circuit, or some other suitable 
device for controlling devices and for signal processing 
of signals. 

Figure 8 shows a signal-conditioning device 8 for the 
10 recorded voltage signal. The measured voltage is signal- 
matched by a resistive voltage division, whereby the 
correct input signal level is obtained for the following 
step which is an anti-alias filter 10. 

15 Figure 8 also shows a signal-conditioning device 9 for 
the recorded current signal. The signal level for the 
current channel is matched to the instrument, either via 
a low-ohm shunt, with the voltage drop across the shunt 
being amplified and becoming the input signal level for 

20 the following step which is an anti-alias filter 11, or 
alternatively, the current signals can be obtained from 
current tongs connected to the instrument. 

The task of the anti-alias filters 10, 11 is to prevent 
25 weighting distortion which arises if the recorded signal 
has a frequency content that exceeds half the sampling 
frequency (see theory for the sampling theorem) . The 
anti-alias filters can be implemented in the form of an 
analogue Sallen-Key low-pass filter according to Figure 7 
30 and have a limit frequency corresponding to half the 
sampling frequency. 
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The level-matched and filtered signals are sampled in 
sampling devices 12, 13 with a sampling frequency, for 
example 6400 Hz. The digital raw data in the form of 
sampled amplitude values is saved in a measurement memory 
5 14 and later constitutes the input data for the 
measurement methods described above. 

Software that controls the signal processor 7 is found in 
a program memory 15. The finished result, that is the 
10 flicker power with sign value and internal impedance, can 
be shown both numerically and graphically in a display 
device 16. The display device can be any known device for 
numerical and graphical display, for example a VDU . 

15 Figure 9 shows schematically a network comprising a 
flickering source Fl, a load LI, a generator G for 
generating an alternating voltage. In Figure 9, the 
current direction for the current I is shown by an arrow 
in the connecting lines. The current goes from the 

20 generator G to the flickering source Fl and the load LI. 
Figure 9 also shows a first measurement point Ml. In 
association with Ml, there is marked a point above 18 the 
measurement point Ml and a point below 19 the measurement 
point Ml. Figure 9 also shows a second measurement point 

25 M2 . In association with M2 there is marked a point above 
18 the measurement point M2 and a point below 19 the 
measurement point M2 . The flickering source Fl emits into 
the network a low-frequency amplitude variation that 
propagates in the direction shown by a solid-line arrow 

30 20 in Figure 9. 
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In the method according to the invention described above, 
the flicker power n indicates in which direction the 
flickering source is located in relation to a measurement 
point. According to an embodiment of the invention, the 
5 sign value of the flicker power is negative when the 
flickering source is located below 19 the measurement 
point and positive when the flickering source is located 
above 18 the measurement point. 

10 At the first measurement point Ml, a negative sign is 
thus obtained for the flicker power n, as the flickering 
source Fl is located below 19 the first measurement 
point. This is because the low-frequency variations in 
current and voltage are in antiphase at the first 

15 measurement point Ml. 

At the second measurement point M2 , on the other hand, a 
positive sign is obtained for the flicker power n, as 
the flickering source Fl is located above 18 the second 
20 measurement point M2 . This is because the low- frequency 
variations in current and voltage are in phase at the 
second measurement point M2 . 

Figure 10 shows schematically a diagram of the flicker 
25 power II for a number of sampling points n during a 
certain period of time when recordings were made of the 
modulated current signals and voltage signals i (n) , u(n). 

Figure 10 shows a first curve Kl on the negative lower 
30 part of the diagram. The diagram also shows a second 
curve K2 on the upper part of the diagram. The first 
curve Kl corresponds to a power signal which, after the 
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creation of the average value, gives rise to a flicker 
power n with a negative value and thus corresponds to 
the flicker power n that was obtained at the first 
measurement point Ml in Figure 9. The second curve K2 
5 corresponds to a power signal which, after the creation 
of the average value, gives rise to a flicker power n 
with a positive value and thus corresponds to the flicker 
power n which was obtained at the second measurement 
point M2 in Figure 9. 

10 

Figure 10 also shows that Kl inverted about the sampling 
axis n corresponds to K2 . This has been demonstrated in 
experiments and the explanation is that the flicker power 
IT changes sign when the propagation of the low-frequency 

15 flicker is changed from going against the power direction 
of the basic tone to going with the power direction of 
the basic tone, or vice versa. It is thus necessary to 
know in which direction the generator and load 
respectively are to be found in relation to the 

20 measurement point in order for it to be possible to 
interpret the sign value in the correct way. The part of 
the diagram that shows the period between zero and where 
Kl and K2 respectively start, shows a period of time when 
the flickering source Fl is not connected in. 

25 
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